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(ABSTRACT:) Electron density measurements have been made in steady-state plasmas in a spherical 
inertial electrostatic confinement (IEC) discharge using microwave interferometry. Plasma cores 
interior to two cathodes, having diameters of 15 and 23 cm, respectively, were probed over a 
transverse range of 10 cm with a spatial resolution of about 1.4 cm for buffer gas pressures from 0.2 
to 6 Pa in argon and deuterium. The transverse profiles are generally flat, in some cases with 
eccentric symmetric minima, and give mean densities of from * 0.4 to 7x 10 10 cm" 3 , the density 
generally increasing with the neutral gas pressure. Numerical solutions of the 1-D Poisson equation 
for IEC plasmas are reviewed and energy distribution functions are identified which give flat 
transverse profiles. These functions are used with the plasma approximation to obtain solutions 
which also give densities consistent with the measurements, and a double potential well solution is 
obtained which has minima qualitatively similar to those observed. Explicit consideration is given to 
the compatibility of the solutions interior and exterior to the cathode, and to grid transparency. 
Deuterium fusion neutron emission rates were also measured and found to be isotropic, to within the 
measurement error, over two simultaneous directions. Anisotropy was observed in residual emissions 
during operation with non-fusing hydrogen- 1. The deuterium rates are consistent with predictions 
from the model. 

[PACS: 52.25.Dg, 52.27.Aj, 52.70.Gw, 52.58.Qv, 52.65.Ff 29.25.Dz] 


I. INTRODUCTION 

The inertial electrostatic confinement (IEC) of 
plasma was originally conceived as a method for producing 
controlled nuclear fusion. It is of research interest in the 
production of high-energy neutrons 1 " 6 and other particles 7 and 
in the study of plasma dynamics, 3, 8-18 as well as in the 
production of energy. 19 " 24 As a neutron source, IEC is of 
interest in detection, imaging and research applications, 
including assaying and tomography, 1, 4 and has been proposed 
for the engineering of sub-critical fission reactors. 5 The rate 
of fusion energy production as a fraction of input power of 
these devices is typically of order 10" 9 , 1 ’ x 4 ’ l9, 2 °’ 25-27 with 
values up to ~1G" 5 reported. 19 In an (ion-injected) IEC 
plasma, the ions travel in radial trajectories and circulate 
through an open spherical, or cylindrical, cathode, and the 
electrons circulate interior to the cathode. The radial 
distribution of density, particularly interior to the cathode, 
which results from these motions is of key significance to the 
fusion and more generally to the plasma dynamics. In some 
cases, there is experimental evidence of spatial oscillations in 
this density, 19 * 28 " 30 and predictions of oscillations, in both the 
density and the electrostatic potential, have played a seminal 
role in IEC research. 

Electron density profiles from D 2 in a pulsed 
cylindrical IEC device were measured by Meeker, et al. 9 
using an intra-cavity laser heterodyne interferometer. 28 
Densities from = 0.3 to 8 x 10 n cm" 3 were measured. 


integrated over a path of 20 cm, and spatial oscillations in 
excess of 30% the mean density were observed. In a previous, 
similar, experiment, Dolan, et a/., 13 using quasi-resonant 
interferometry of metastable neon atoms, determined a 
ceiling of 2xl0 ,! cm" 3 for the election density and estimated 
specific values for this density based on a theoretical 
dependence of the metastable atom density on the electron 
density. Dolan, et al. y also presented calculations of intra- 
cathode density and potential profiles using a 1-D kinetic 
model and found solutions predicting densities consistent 
with the neon experiments. These particular solutions used 
step function distributions of total energy and angular 
momentum, and the same basic model has been examined by 
other workers, 14 * n * 29 ’ 31-33 using also delta-function 
distributions 14, ,7 * ,9 * 33 and quasi-thermal distributions, 17 * 32 
typically for the purpose of testing assumptions about the 
physics or for comparison with experiment 

We report here electron density measurements made 
in spherical steady-state IEC Ar and Dj plasmas using 
microwave interferometiy. The line-integrated, or areal, 
densities are from < 1 to s 8 x 1G 11 cm -2 , over paths of 15 
and 23 cm, corresponding to uniform densities from = 0.4 to 
7 x 10 10 cm" 3 . Transverse scans of the plasma over a distance 
of 10 cm were performed, and the resulting profiles are 
generally flat with oscillatory structure evident in some of 
the data. The spatial resolution for these measurements was 
limited by the beam diameter to about 1.4 cm. By 
comparison, the laser heterodyne technique has a spatial 
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resolution of « 1 mm, but is not amenable to steady-state 
discharges, and is perhaps less sensitive due to die short 
wavelength. In Sect. IV, calculations from the 1-D 2-species 
kinetic model are summarized and specific theoretical 
transverse profiles are compared with the measurements. A 
modified rectangular distribution function and the thermal 
equilibrium distribution are introduced, largely to produce 
flat central profiles, and explicit connection to the 
interelectrode region is made. 

Measurements of the D-D fusion neutron production 
rate have also been made using a pair of calibrated 3 He 
detectors. The measured rates are consistent with rates 
reported in the literature 4, 25-27 and are compared with 
predictions from the model. The dual-channel configuration 
is used to quantify the neutron emission isotropy over two 
directions, and comparative studies using both H 2 and D 2 
provide information on the cathode fusion contribution, 4 
which is less isotropic. Both the electron density and the 
neutron emission measurements were made using two 
cathode diameters for a fixed anode. 


II. EXPERIMENT 

The IEC apparatus is comprised of a 60-cm diameter 
spherical vacuum chamber with an interior concentric 
cathode grid and supporting hardware. The vacuum chamber 
wall is grounded and forms the anode of the circuit, and 
voltage is applied with a 100-kV, 5 kW power supply. The 
voltage V g and current Is are measured independently of the 
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FIG. 1 . Schematic of IEC device showing (a) basic configuration and (* add 
V and I meas ?; H20 in/out?) (b) microwave and He3 diagnostics. 
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FIG. 2 IEC operating voltages as functions of the Paschen product pd for 
protium and deuterium using two cathodes. 


supply using an external voltage divider and an integrated 
parametric current transformer, respectively. The flow of 
neutral gas through the chamber is controlled by a 0.1 slm 
MKS flow controller, and the chamber pressure p is set by the 
flow and pumping speeds. The chamber is double-walled and 
water is circulated through the inter-wall cavity for thermal 
regulation. This regulation was found to significantly enhance 
the stability of operation, particularly at high voltages. A 
schematic of the IEC device is given in Fig. la. 

The cathode grid consists of 7 “great circle” stainless 
steel loops, which are pre-fabricated o-ring seals, welded 
together to form a symmetric lattice of triangles of two, 
similar, sizes. Two diameters, 15.2 and 22.9 cm, of cathodes 
of this design were used. In each case the wire loop cross 
sectional diameter is 0.0794 cm, resulting in geometric 
transparencies of 0.964 and 0.976, respectively, for the two 
cathode sizes. A second 15 cm grid design, based on 
“latitude” circles, was also used for some of the 
measurements. This “L-grid” consists of three great circle 
loops with mutually orthogonal axes, of which two loops are 
each flanked by symmetric pairs of smaller, coaxial, circles; 
the geometric transparency is = 0.95 (*). 

The experiment was operated in discharge mode 20 so 
that the operating voltage is fixed by the neutral gas chamber 
pressure. Once breakdown has occurred, small changes in the 
voltage suffice to adjust the current over the full range of the 
power supply. 3 Argon and deuterium gases were used for the 
interferometry. Deuterium was used for the fusion work and 
hydrogen-1, or protium, was used for comparison with 
deuterium. The experimental conditions for these gases are in 
the large mean-free-paih regime of the Paschen curve for 
breakdown voltage. 34 Fig. 2 gives Paschen data for the two 
cathode sizes for protium and deuterium, taken from 
experiments measuring neutron emission rates. The striking 
distance d is the interelectrode spacing, 3 and it is seen that the 
two cathodes fall on different curves. 

The neutral gas particle density n 2 is required for the 
theoretical fusion reaction rates in Sect. V, as well as for 
calculation of the ionization fractions (Sect. VI), and the n 2 
used here are taken from the pressure p assuming a nominal 
neutral temperature T 2 . This temperature T 2 is generally 
higher than the ambient temperature T** and depends on the 


2 







fraction a of the electrical power | V g 1^ channeled into the net 
heating of the neutral gas. In particular, the change in 
temperature experienced by the neutral gas as it passes 
through the vacuum chamber is 


AT |F < 17 * 

A T, =a—£ 

2 Mpk 


( 1 ) 


where M is the rate at which neutral particles enter (or 
leave) the chamber, known from the flow controller, is a 
constant in the range of from 5/2 to 9/2, depending on the gas 
species and temperature, and k is the Boltzmann constant. 
The fraction a is not known, but unless it is very small, 
substantial heating will result; the value required for A T 2 = 
b is a £ 0.002 for typical microwave data in Ar and a = 
4x1 O' 5 for typical fusion data. 


111. MICROWAVE INTERFEROMETRY 

A three-frequency microwave interferometer was 
used to make electron density measurements in the IEC 
plasmas. Optical access is provided by a pair of opposed 
viewports in the equatorial plane of the vacuum chamber. 
These viewports have a clear aperture of 10 cm and the 
separation of the windows along the optical axis is 
approximately 76 cm. The windows were made of either 
alkali borosilicate (Coming 7056) glass or plexiglass. The 
transverse spatial distribution of the microwave beam 
irradiance is approximately Gaussian with a nominal 1.25 cm 
focal spot. The two microwave frequencies used for the IEC 
work were 70 and 90 GHz, and the full-width at half- 
maximum of the 70 GHz beam waist was measured to be 1.8 
cm. The transverse scanning range for the microwave beam 
was 10 to 12 cm. Stepper motors were used both for the 
transverse scanning and for the phase calibration, which was 
performed by direct axial translation of the interferometer 
transmitter. A schematic of the microwave configuration is 
shown in Fig. lb and further details are given in Ref. 35. 

The interferometer was calibrated at each location in 
a transverse scan. Just prior to, or immediately following, the 
plasma measurement, the transmitter was translated axially 
over a range of about one fringe centered on the measurement 
position, which was chosen near a zero-crossing. The 
interferometer has quadrature capability and the data 
acquisition was performed off-quadrature, with shifts of < 
55°, so that there were two useful channels of data at each 
transverse location, both of which were used in the data 
reduction. Also, in many cases, calibration scans were 
performed with the plasma present in the chamber for 
comparison and as a check. The axial scans were performed 
in discrete steps and the position of the transmitter was 
recorded at each step, and this resulted in typical calibration 
times of = 100 sec due to constraints on the polling speed for 
the stepper motor (the scan itself requires only a few 
seconds). The plasma measurements themselves were 
performed by cycling the plasma on and off several times, 
with each cycle having a duration of approximately 10 
seconds, and the phase shifts were subsequently determined 
by averaging the data flanking the termination points in these 


cycles. The operation of the IEC discharge is stable over 
periods of many minutes, particularly at low voltages, and the 
restriction of die measurement time in this way was done 
largely to reduce the measurement error from drift in the 
interferometer signals. 

A Monte Carlo error analysis was performed for the 
phase shift measurements, based on the individual calibration 
files and the measurement data itself. 35 The resulting error bar 
for the measurement of electron density, as presented here, 
represents an estimate of the unit width ( V2 standard 
deviations) of a Gaussian distribution. This density error is 
generally substantial, as a fraction of the measurement value, 
due to the small phase shifts, which are typically of order 10~ 2 
fringe. 

Unlike the cylindrical case, the interferometer beam 
in spherical geometry necessarily has a radial component and 
the beam must pass through the cathode grid. In the present 
case, diffraction and absorption due to the grid are expected 
to be minor. The grid apertures are several centimeters in 
dimension, and the ratio of the diameter of the wire loops to 
the interferometer wavelength is < 024. In this connection it 
is of note that the microwave polarization was perpendicular 
to the transverse direction, so that any wire loop with an 
approximately transverse orientation, where it might have 
maximum effect, would tend to select for the beam 
polarization. In any case, any effect of the grid on the beam is 
present during calibration as well as during the 
measurements; also, trial data reductions using the 
calibrations made with the plasma present gave very similar 
results to those made in the absence of the plasma (apart from 
the sign of the phase shift). 

The spatial resolution for the interferometry, set by 
the beam geometry, is taken to be 1.8 cm for the 70 GHz 
data, from the waist measurement, and 1 .4 cm for the 90 GHz 
cases, from wavelength scaling and the manufacturer's 
specification. The wavelengths for these two beams are 0.429 
and 0333 cm, respectively, and thus there are transverse 
spatial gradients in the microwave irradiance on scales 
comparable to the wavelength. Moreover, there are 
potentially gradients in the electron density on scales 
comparable to the wavelength, due to intra-cathode radial 
structure. The effect of such density gradients on the 
microwave detector signals is generally difficult to calculate, 
particularly given the irradiance gradients, and in the data 
reduction here this effect is neglected based on the 
assumption that any error thus introduced is small in 
comparison to the phase shift error. The inferred phase shift is 
taken as proportional in the usual way to the line-integral of 
the density, i.e., to the product of the path length and the 
mean density in the volume occupied by the beam. 

For the interferometry, the apparatus was operated 
with cathode voltages = 1 to 5 k V and power supply currents 
= 50 mA; neutral gas pressures were from = 02 to 6 Pa. 
Visually, the IEC plasma under these conditions typically 
manifests a fairly well defined core region with diameter 
comparable to, but less than, that of the cathode. Outside the 
cathode, the electron density is necessarily very low; e.g., for 
50 mA of direct electron flow from the 15 cm cathode to 
anode at 1 kV, the mean density is s 10 5 cm -3 . Therefore, the 
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plasma may be considered as confined to the interior of the 
cathode. 

Transverse profiles of the areal electron density N 
from the Ar plasmas are given in Fig. 3. These data were 
collected on several days over a period of a few months. The 
profiles are generally flat across the center of the plasma and 
have average densities of from = 4 to 8 xlO 11 cm -2 . Most of 
the radial variation is comparable in magnitude to the 
measurement error, but there is a significant tendency toward 
eccentric symmetric minima in some cases. This is clearest in 
the standard 15 cm grid data, for which distinct minima are 
evident at x s ± 4 cm. It should be noted that clipping of the 
microwave beam by the viewport aperture is present in some 
degree for most of the data. At ±5 cm, where the beam is 
nominally bisected by the window edge, attenuation of the 
signal is > 80%, and this is reflected in the error bars. The 
transverse step size was 0.635 cm, except for the L-grid data, 
most of which was taken with steps half this size. The 70 


x{cm) 

-6 -4 -2 0 2 4 6 



FIG. 3. Ar interferometer transverse scan data. The left axis gives areal 
electron density N in units of 10 12 cm* 2 , the right axis, phase shifts in 
degrees. 


GHz microwave frequency is somewhat more sensitive to the 
density than the 90 GHz component, but die accuracy of the 
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FIG. 4. D 2 interferometer transverse scan data. The left axis gives areal 
electron density N in units of 10 12 cm' 2 , the right axis, phase shifts in 
degrees. 
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FIG. 5 Pressure scan data for the interferometer. The left axis gives areal 
electron density N in units of 10 11 cm' 2 , the right axis, phase shifts in 
degrees. 
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measurements is affected also by the signal strength and 
calibration quality. The transverse scan data for D 2 , collected 
over a comparable time period, are presented in Fig. 4. These 
profiles are qualitatively similar to those in Ar, but the 
densities are distinctly smaller, in the range of from = 0.7 to 3 
xlO 11 cm -2 . 

Pressure scans were also performed at the center 
location (x = 0). This data is given in Fig. 5. The electron 
density generally increases with neutral gas pressure, but 
there are indications of limiting densities at the higher 
pressures. This is clearest in the Ar data for the 23 cm grid, 
which shows a constant density in the range 4-6 Pa. This 
behavior is reminiscent of that of fusion rate density at lower 
pressures observed by Hirsch and Miley, 20 in which the 
fusion rate increases with increasing pressure only up to a 
point (about 1.3 Pa), and is discussed in Sect VI. 

For the purposes of comparison with the model in 
the following section, the essential features of the 
experimental transverse profiles are that (1) the areal 
densities are of order 10 u to 1G 12 cm -2 , (2) the densities are 
approximately constant in the center, giving flat profiles, and 
(3) some of the 15-cm grid data show minima in the densities 
at * = ± 4 cm. 


IV. STEADY-STATE 1-D KINETIC MODEL 


The convention adopted here for the definition of 
space charge currents includes charges traveling in both 
radial directions, inward and outward, 13, 19 rather than in a 
single direction only; 14, 36 this total current corresponds 
directly to the number density of charge carriers, which 
drives the Poisson equation. In particular, the circulating ion 
current 7 0 is related to the ion current collected by the 
cathode I, r by 37 


l ic 

_Al = 2 r e 

he i -r\ 


( 2 ) 


where y e is the effective cathode grid transparency. y e is 
related to the open fraction y 0 = l-cr 0 1 A , where cr 0 is the 
cross-sectional area of the grid and A is the surface area of 
the grid sphere, by the grid absorption efficiency. 


n = 


i-ro 


(3) 


i.e., rj = a g lcr 0 , where a g is the cross section for 
absorption of ions by the grid. Dolan has given an expression 
for y € /y 0 which depends on the grid wire spacing and 
diameter and for which a typical value of 77 = 1.35. 37 On the 
other hand, Miley, et a!., have given fits of empirical data to 
Paschen curves, supported by theoretical arguments, which 
give typical 77 ^ 0.1. 3 

The current which is measured experimentally is the 
power supply current I s , and it is related to the circulating 
current by 


4> 


i+s 1 - Y y 


(4) 


where 5 is the secondary electron emission coefficient and 
I hv is the secondary photoelectric emission rate. For the 
nominal operating current I s = 50 mA, values of / 0 , along 


Tabic I. Geometric transparencies and nominal circulating currents for the 


two cathode grids (/$ = 

50 mA, 17 * 1). 
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/o 

15 cm grid 

0.964 

1.24 A 

23 cm grid 

0.976 

1.87 A 


with y 0 , are given in Table I for the two cathode grids for the 
case 77 = 1, using, from the interferometry voltages of 1 - 2 
kV, S = 0.1, 28, 38 and assuming 1^ « I s . These values of 
I Q correspond also to 77 = 0.56 for the case 5 = l. 39 Also of 
interest is the electron space charge current into the anode 
I e ; this current flows in one direction only and, although I s 
has other contributions, the net result is 


h = h- (5) 

It is these primaiy currents J 0 and I e , then, which 
anchor the steady state 1-D model. This model consists of the 
Poisson equation, 


-2 


—{ 
dp V 


,2 df 

dp j 


= ”o-”l 


( 6 ) 


and a prescription for the distribution over energy of the 
charge carriers. 19, 12 ‘ 14, 17, 3133 Here p and ^ are the 
dimensionless radius and potential, respectively, defined as 


P = 


r s 


4 = JL 

* v* 


(7) 


where r and V are the radial position and electrostatic 
potential, respectively, and the subscript g denotes values at 
the cathode grid; /Iq and n x are the number densities of 
charge carriers from the ^ = 0 and 1 electrodes respectively; 
and 



a 

Ie 


5.53x10* 


(r g cm -1 ) 2 


cm 


( 8 ) 


is a characteristic density, in which is die free-space 


permittivity and q e the electronic charge. It is of note here 

that this formulation is applicable to grid voltages of either 
polarity, 13 and that it is convenient in some calculations to 
reference the potential to the anode. 19 However, for the 
explicit treatment here: V g <0; = 0 at the anode; the 

subscript “0” will denote ions; the subscript “1” will denote 
electrons; and the index i = 0,1 will be used to denote either 


(or both). The ion charge is taken to be q 0 = -q x =\q e \. We 

introduce two dimensionless currents, the mono-energetic 
perveance, 19 


K, = 


*x^\Vj n \\2q e \ 


m. 


xl/2 


(9) 


and the general form 



( 10 ) 


where C, is a distribution function normalization factor, 13 
and current-normalized densities n, , which are related to the 
Poisson densities by 


hi = k; x ^- 


(ii) 
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The procedure then is to form a pair of functions 
and n x (p,$) , considering p and ^ as independent 
variables, by specifying constraints on the distribution over 
energy of each species of particle. This fundamentally 
provides a relation between <j> and the equations of motion. 40 
Defining the (conserved) angular momenta, 

A = (12) 

the (conserved) total energies, 

£, = 9/ V(r) +^-m i ulj + ^ 2 , (13) 

2. 2m i r 


the normalized energies and momenta, 

« - -a 

<ii v g 


, (14) 

2m, r g E g 

and the integral-normalized distribution functions / 0 and /, , 
the normalized densities are 


n,(p,4) = fafaf,(j>,4,e n X,), (15) 


where: u r and u t are the radial and transverse velocities, 
respectively; E g = \q e V g \; and the e , , ,4 are to be implicitly 
constrained so that the Kinetic energies, #-e 0 and e x 
and the radial kinetic energies \$~-e i \ -A?p~ 2 , remain 
nonnegative. 13 Using these functions n, , solutions #(p) of 
the Poisson equation are obtained for values assigned to the 
A, resulting in density profiles n,(r) = n g k t n,[p,^(A,/>)] . 
In the context of kinetic theory, the distribution functions f 
provide solutions to the Vlasov equation. 17 

It should be noted that while the solutions of the 
Poisson equation depend on the circulating current through 
Eq. (9), these solutions may be obtained without explicit 
specification of the current, cm- of the cathode voltage or 
radius. Given a solution ^(p) , the physical number densities 
n t require specification of V g and r g , but not m , , and the 
circulating current / 0 requires specification of V g and , 
but not r g . The relation between / 0 and n^(r g ) is 

A 

4 * r gM , ( 16 ) 

n g ^0 fao fa ^ O )u r o(l,l,£ : 0’'^o) 

and that between I e and w,(r 0 ) 

u 

4 / rr ° h \ , (17) 

n g k x fa x f x (aA£\,^)v rX {aA£\,*\) 

where r a is the anode radius and a = r a /r g ; the 
normalization facto rs C, are evaluated from these 

expressions as C, = <j2E g /m t f v r4 . 13 

In the simplest case, the particles are mono-energetic 
and have purely radial motion. This is the original model due 
to Hirsch and Porter. 12, 19 All ions have £ 0 = 0, all elections 
havef, = 1 , and the densities are simply related to the space 
potential through conservation of energy, becoming singular 
at the source electrodes, where the velocities vanish, and at 
the origin. In this case, there are clearly defined limits, 
imposed by the space charge, on the magnitude of the 
currents. Given a ratio of electrode radii a, there is a unique 
solution for which both currents are space-charge limited. 8, 19 



FIG. 6. Mono-energetic solutions of the IEC Poisson equation. 


For the 15 cm grid, a = 4 and this limits the ion current 
between r a and r g to K 0 = 0.296, corresponding to 2.29 mA 
for Ar ions at V g = -1 kV. For the experimental case, 
however, I e = 50 mA ( K x = 0.024), giving K QumL - 0.0949 
and 7 0 = 0.734 mA. In either case, K 0 is continuous across 
the cathode, and K x , which need not be, takes up a new value 
interior to the cathode based on this AT 0 , and so on. 19 The 
densities and potential for the case of space charge limit on 
both currents, and for the case of limiting currents subject to 
l t = 50 mA, are given in Fig 6. The free-space density 
singularities interior to the cathode have been called virtual 
electrodes, and it is of note dial, although the /?, become 
singular at all source electrodes, these singularities are 
integrable. Hie mean density in a spherical shell 
\P~Ps l« A » where p s is the radius of the singularity, is 

n,(p-p,)= , , 4 f u £ dp' p' 2 n,(p') , (18) 


in which /i,(p) = n K t p 2 & 1,2 and 
s (9K t / 4) ' (I p - p, I / p, ) 4/3 ,*• 19 giving 


"Ap-p,) 


•JUK, 

p\\p-p. 



(19) 


Eq. (19) is used in computing the mean densities for Fig. 6, 
given in Table II below. Also of note is that the data plotted 
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in the figure are cut off arbitrarily at small p ; in the model, 
successive interior singularities exist for indefinitely 
decreasing p (and increasing ratios of virtual electrode 
radii). 


The mono-energetic space charge limit on density 
may be effectively circumvented by dispersion of the particle 
energy, and these dispersed distributions also generally 
remove the electrode singularities. Rectangular distribution 
functions, based on the Heaviside unit step function h y have 
been examined by Dolan, et a/., 13 and others. 14, 17, 31 ~ 33 The 
case 13 

/*</>.#)« 

— ? Tl | <y _}-|-l/2> (20) 

n g ki /*(*„ Dp 2 \\<Kp)-e\-tfp ] 

where 


/?,(*;, 4 ) s 

h(€i -e iM max "4) ’ 

is plotted in Figs. 7(a) and (c) for the 15 cm grid, using 

^O.min ^l,m« ^ 


( 21 ) 

( 22 ) 


Physically, Eq. (20) represents an even distribution of ion 
birth sites. The particles take up orbits in which they oscillate 
between radii which share a common potential. Thus the 
range of particle energies at a point is constrained only by the 
local potential, i.e., low energy ions are absent where 1 - ^ is 
large and low energy electrons are absent where <f> is large. 
However, in contrast to density profiles from these 
distribution functions given for cylindrical geometry in Ref. 
13, the n ( from Eq.(20) in Fig. 7 remain singular at the 
origin, increasing like p~ l (^ constant). This singularity 
may be removed by thermalizing the angular momentum. 
Mathematically the singularity is due to convergence of 
particles at the origin and these particles must have 4 -> 0. 
If the distribution function vanishes for ^ -► 0, which is a 
property of thermal distributions, the n t (0) remain finite. 
This is illustrated in Fig. 7 as well, which also plots solutions 
using the parabolic distribution function, 

(23) 

where C 0Ji is given explicitly in Ref. 13 (there called C,), 
and C, H and C i L are derived in the same way [from Eqs. 
(15) and (16)], TTie step function limits are those of Eq. (22). 
The distribution over kinetic energy for Eq. (23) is sensibly 
identical to that of thermal equilibrium at high temperatures, 
except that it is truncated by the step functions where the 
kinetic energy equals the local potential (1-^ for f XyL ). The 
mathematical form of f i4 X e i) - jXi(^,4)^4 is, within 
the limits, generally parabolic, although for both the 
rectangle and the parabola functions, limits on 4 may 
modify the shape. 

In Fig. 7, as in the experimental case of Fig. 6, the 
currents have been maximized subject to the constraint I s = 
50 mA. For the dispersed distributions, however, the limiting 
currents are not defined as precisely as with the mono- 
energetic case. The general procedure adopted here to 
determine a value for k iaax is to increase k t until the 


Kan) 

II 15 19 23 27 30 



1<T 

10 # 

10 6 | 

10 4 

10 *° 

10 8 

10 6 I 
10 4 

1.0 

0.8 

0.6 ^ 
0.4- 
0.2 
0.0 


FIG. 7. 15 cm grid density and potential profiles for kinetic solutions using 
dispersed distributions (rj = 1, £-0.1). (a) gives for the rectagular 
functions Eq. (20) and the parabolic functions Eq. (23); and (b) gives n, for 
the thermal functions Eqs. (24) and (25). The two cases in (b) are for the 
superposition of broad (full) and narrow (clipped) step functions on die 
thermal distributions, (c) gives 1 for all cases. In all cases, / s = 50 mA, 
V g = -1 kV, and the electron and ion densities are essentially equal for 

p <= 0.8 . 


numerical integration of the Poisson equation begins to foil, 
where a failure is defined in terms of a merit function for the 
integration. This merit function is based on the degree to 
which the target boundary conditions are met and die degree 
to which the target boundaries (in p) themselves are 
attained, and an arbitrary upper limit on the merit function of 
= 0.01 was used. The numerical integration is performed by 
the shooting method using an adaptive step size Runge-Kutta 
algorithm, 41 and the integration proceeds either in a single 
radial direction, starting from an electrode, or, for high 
densities, in both radial directions, starting from a variable 
inter-electrode or intra-cathode p . The boundary conditions 
are ^'(0) = <f>(a) = 0, #1) = 1, where the prime denotes 
differentiation with respect to p . In Fig. 7, then, Eq. (5) is 
used to set k x for p > 1 , Table I is used to set (same for 
all p ), and iterative integration is used to maximize k, for p 
< 1. The resulting current parameters and mean densities are 
given in Table 11, along with values from analogous 
calculations for the 23 cm grid. To the extent that the exterior 
region is not of interest, the only requirement imposed by the 
exterior k, is that it be compatible with k 0 for a given a. 

For comparison. Figs. 7 (b) and (c) give profiles for 
the thermal distribution functions: 

yj r (£,,/$,) = 


JLn'hXeAHp 1 
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where 



(24) 


(25) 


and 7V and T (J are the radial and transverse temperatures, 
respectively. In the figures, and in all the work here, 
£ 0 = T ( 0 and T r , = T t x are set to give mean kinetic energies 
LJ i , at a given p^ y equal to those of the rectangular 
distributions Eqs. (20-22): 

~ • (26) 
Two version of the thermal functions are plotted, one using 
“fiilP thermal distributions with e Qjniak = -1, e l amx = 2, and 
one using truncated, or “clipped,” distributions with e Q nua = 
0, e l mMx - 1 ; in both cases, the limits are otherwise as given 
in Eq. (22). For the full thermal distributions the limits are 
sufficient, for die temperatures here, that U i = 3kTJ2 to 
good approximation, and for the clipped distributions, a 
bisection routine is used to find the temperatures, generally 
giving T t > 2 U i /(3k) . The currents and mean densities are 
listed in Table II. 

Physically, the full thermal distributions have 
particles with energies unbounded by the local potential, but 
this is largely inessential to the treatment here. The fraction of 
the ion population with U 0 > \ q € V | , or of electrons with 
U x > E g -\q e V\, is typically « 0.1 and, as shown in the 
figure, die results are very similar for both the clipped and 
full distributions. The key feature of the rectangular n t (p f ^) 
is that they decrease as ^ approaches the potential of the 
source electrode ( n 0 decreases away from the cathode, and 
vice versa) because the range of kinetic energies available 
decreases. The decrease in kinetic energy range occurs at 
high energies and consequently U t also decreases. In 
preserving U f , this feature of the rectangular distributions is 
transferred to the thermal distributions, and as a result, the ion 
temperature decreases and the electron temperature increases 
in moving away from the cathode (potential). 

The densities for the distributions in Fig. 7, while 
larger than the mono-eneigetic densities, remain insufficient 
to explain the microwave results. Transverse profiles from the 
solutions of Table II are given in Fig. 8. These profiles are 
based on a model of the microwave beam having a step- 
function radial profile and a Gaussian-beam axial profile; the 
central radius, or “waist,” is 0.7 cm and the radius at the grid 
is 1.0 cm, giving a Rayleigh range = r g . The mean areal 
densities from these transverse profiles are of order 10 10 cm' 2 . 
Incrementally higher densities are possible with k 0 = it, ; 
e.g., for the parabolic functions, the maximum current is near 
k Q = k x = 1000 and gives, for the 15 cm grid, n 0 = 2.55x1 0 9 
cm' 3 , n x = 2.03x1 0 9 cm' 3 (cf. Table II), I Q = 1.93 A , r\ = 
0.357, and the external solution is valid for l s - 50 mA. 
However, the ion and electron densities from these kinetic 
solutions are systematically disparate where the electric field 
is strong, near the grid, and, as pointed out in Ref. 28, such 
disparity is limited by the cathode voltage. In particular. 
Gauss’s law requires (*) 


Table n. Interior parameters (p< 1) for Figs. 6 and 7, and for calculations 
for the 23 cm grid analogous to the cases in Fig. 7. The mean densities n, 
are in units of 10 6 cm' 3 . For the mono-energetic case, C# = 1, and I o = 



para- 

mono- 

rectangular parabolic 

clipped 

thermal 

Grid 

meter 

energetic 



thermal 



*0 

0.0949 

241 

641 

1544 

2303 

15cm 

*1 

0.00987 

65 

1300 

2300 

2750 



245 

1529 

2140 

2140 

1682 



245 

1051 

1770 

1729 

1276 


*0 

0.250, 

0.00565 

363 

967 

2329 

3473 

23cm ~ 


0 0567 

20 

850 

1500 

3100 



3.97 

462 

1044 

1014 

984 



3.97 

175 

813 

761 

753 


J [dpp - 2 £ dp' p ' 1 [/*,(//)-/»,(/>')] S . (27) 

For the 15 cm grid at 1.74 kV, for example, the uniform net 
density limit is = 10* cm" 3 . Also of note in Fig. 8 is that for 
all of the kinetic solutions N decreases monotonically with 
increasing x. 

To further increase the densities, then, we explicitly 
invoke, to varying degrees, the plasma approximation, 
replacing the kinetic electron density n, with a quasi-neutral 
density h x defined as 

h x - (l-K)no + vni, (28) 

where v is an arbitrary neutrality parameter. For v 0 and 
k 0 =k ]y Eq. (28) reduces to the conventional plasma 
approximation. 40 However, while using n, with the 
distribution functions Eqs. (22)~(24) permits high intra- 
cathode densities, continuity of / 0 at the cathode requires 
also high ion density exterior to the cathode, where the 


x(cm) 


0 2 4 6 8 



FIG. 8. Transverse profiles of kinetic solution data for (a) the 15 cm grid 
and (b) the 23 cm grid. All profiles except the mono-energetic solutions are 
referred to the left ordinate axis. The units are 10 10 cm -2 . 
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plasma approximation is untenable because of the constraint 
on I e . These high exterior densities are ordinarily associated 
with small values of a , 8 and this ratio is fixed here by the 
experiment. Mathematically consistent solutions may be 
formed in this case by £ 0min > 0, which, physically, removes 
those ions with sufficient energy to approach the anode and 
increases the concentration of ion space charge near the 
cathode. 

For the thermal functions, the high currents at the 
cathode may be markedly reduced while preserving the 
interior densities by truncation of the distribution at die low 
energy limit, but this truncation is much less effective for the 
parabolic functions. This is significant because, for both the 
thermal and the parabolic functions, high currents at the 
cathode result in small values of 77, from Eq. (3). The 
truncation e Q max < 1 is thus used here to calculate low-current 
thermal solutions with 77 * 1, a relation discussed further in 
Sect. VI. 

In treating the quasi-neutral cases, then, the energy 
and angular momentum limits will by default be those of Eq. 
(22) except that e Q mm = 10~ 3 , and, for the thermal case, 
£ l mx - 2. Where the limits differ from these values, as with 
the low-current thermal solutions, they will be given 
explicitly. Also, for the exterior solution, because of the high 
ion densities some measure of neutrality js imposed on the 
thermal electron currents by using U l T = U Q H in n i r (p>l). 

Fig. 9 gives radial and transverse profiles from 

r(cm) 


012345678 



FIG. 9. Radial (a) and transverse (b) profiles for the 15 cm grid in the 
plasma approximation. V g = -1.48 kV, v= 0.004 in all cases. £ 0mmx - 0.65 
for the low-current thermal case. Shown in (b) only is a second low-current 
thermal case (broken line) with 7 = 1 ( £ 0 toax = 0.6671). The experimental 
case is that of Fig. 3 (b). 


Table III. Parameters for the theoretical cases in Fig. 9. 



thermal 

low-/ thermal 

parabolic 

IHEtfl 

3.223 

1.961 

0.8537 

MM3MSWM 

3.195 

1.960 

0.8500 

/o(A) 

46.7 

0.780 

29.7 

7 

0.027 

1.57 

0.042 


solutions in the plasma approximation and Table IQ lists 
parameters used in these solutions. The radial densities have 
maxima at or near the cathode, and the currents are chosen to 
give mean areal densities comparable to those of an 
experimental transverse profile from Fig. 3. This 
experimental profile is also shown in Fig. 9, where it has 
been folded about a best-fit symmetiy axis near the nominal x 
= 0 and averaged over ±x using cubic spline interpolation. 
Additional folded profiles from the measurements are given 
in Fig. 10, along with other selected solutions from the model 
in the plasma approximation; theoretical parameters are given 
in Table IV. In these cases, the ratio ko / k x >1 and the 
parameter v > 0 are used to determine the shape of the 
transverse profile near the cathode. Although the 
measurement uncertainties qualify the agreement, these 
calculations demonstrate quantitative consistency between the 
1-D model and the microwave measurements, particularly 
with respect to features (1) and (2) listed at the close of Sect. 
3. 

With respect to the eccentric minima, item (3), these 
dispersed solutions do not exhibit local minima, and this is 


x(cm) 

012345678 



FIG. 10 Comparison of experimental and theoretical transverse profiles in 
the plasma approximation for (a) the 15 cm grid and (b) the 23 cm grid, 
^o.nn = 0-65 for the low-current thermal case. The measurement cases, 
identified by pressure, are, from top to bottom, those of Figs. 3(c), 4(a) 
(explicit error bars), 3(d), and 4(b) (explicit error bars), respectively, (label 
error bars) 


Table IV. Parameters for the theoretical cases in Fig. 10. 


P< 1 


15 cm 

23 cm 


low-/ 

parabolic 

thermal 

parabolic 

*□ no 4 ) 

4.302 

0.375 

9.548 

1.4572 

MIO 4 ) 

4.300 

0.360 

9.530 

1.4470 

/o(A) 

1.01 

23.6 

78.0 

37.9 

r? 

1.22 

0.053 

0.024 

0.050 

V 

0.002 

0.01 

0.0002 

0.004 

% ( 10 lo crrf ! ) 

4.064 

1.560 

4.039 

1.822 

wi( 10 10 cm -3 ) 

4.058 

1.506 

4.031 

1.810 
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emphasized by the example in Fig. 9. In fact, the implication 
of the pronounced minima in Figs. 3 (a) and (b) is that there 
are substantial local gradients in the radial density, suggestive 
of a virtual electrode. For this reason, solutions with 
inflections in the potential, like the mono-energetic cases, are 
of interest. However, unlike the mono-energetic cases, these 
solutions must also have flat central profiles and high interior 
densities. The flat center may be retained in a quasi-mono- 
energetic solution by appropriate choice of the limiting 
angular momentum, and the density may then be increased by 
superposing quasi-neutral populations from Eq (28). 
Examples of this are given in Fig. 11 and Table V. The 
potential profiles in the figure are determined by narrow 
parabolic distributions having Ae= 1(T 5 and = 0.3, and 
the densities are largely determined by the thermal 
distributions. For simplicity, ko / k x =1 and v = 0 for the 
thermal functions. The potentials ^ in Fig. 11 are double- 
welled and have local minima at the “virtual anodes” of 4.04 
and 5.701 xl(T 3 for the 1.06 kV and 1.48 kV cases, 
respectively. Thus the virtual anodes are not truly anodes, in 
the sense that the ions which form them are not bound by 
them. The ions have sufficient energy to breach the potential 
barrier and therefore the ion currents must be continuous 
across the anode. For p > 1, the parabolic ion distribution is 
shifted to accommodate a so that, in this case, f 0tL is not 
strictly continuous at the cathode, although f 0 j and I 0 are. 
Also shown in Fig 11 are the folded profiles from Figs. 3(a) 
and (b). The double-well solutions are qualitatively similar 
and have appropriate mean densities, but the experimental 
minima are more pronounced than those of the calculations. It 
is of note in this regard that calculations of the transverse 
profiles of these double-well solutions which use high spatial 
resolution (<1 mm) do not differ dramatically from those 
given in the figure. 

One feature of the dispersed distributions which is 
counter to that found in the mono-energetic case is die 
concentration of ion space charge at the cathode, and this is 
true of both the kinetic and the quasi-neutral solutions. For 
the low-current thermal cases, this peak in the ion density is 
moved radially inward from the cathode, as in the total 
density of Fig. 11(b). It is worth mentioning in this 
connection that calculations from a Berkeley Particle-In- 
Cell/Monte Carlo code have shown also an intra-cathode ion 
density peak, which is attributed to thermal effects, 1 although 
these calculations do not give flat centers (nor exterior 
peaks). Also of note is the intra-cathode maximum in 
cylindrical geometry illustrated in Ref. 13 (cf. Fig. 5); in this 
case, the shift in the maximum, as compared with Fig. 7, is 
due to limits on the angular momentum. More generally, the 
dispersed distributions have the property, with limits on e, 
and 4 tending to exception, that the shape of the curve n 0 (p) 
follows qualitatively that of <f(p\ and h x (p) that of l-^(p). 

V. NEUTRON MEASUREMENTS 

D-D fusion neutron emission rates S were 
measured using two Eberline Swendi-Il 3 He neutron 
detectors. The calibration constant furnished by the 


Kcm) 
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FIG. 1 1 Potential (a), radial density (b), and transverse density (c) profiles 
from double potential well solutions using quasi-mono-energetic parabolic 
functions and low-current neutral thermal populations. v = 0forp<l(v = 
1 for p > 1), Asx = ~ 10' 5 , Ax^mx = 0.3. if - 1 by choice of 

SifTjmx The bands in (c) are measured transverse profiles in Ar from Figs. 
3(a) and (b), and V g and I, for the solutions are taken from these cases. 


Table V. Parameters for the theoretical cases in Fig. 1 1 . 



1 .06 kV solution 

1.48 kV solution 


p< i 

P >\ 

P< 1 

P> 1 

for (10 s ) 

1.86 

0.039xl0" s 

0.620 

0.023xl0' s 

*.l(10 5 ) 

181 

- 

215 

- 

fill min 

0 

0.001 

0 

0.001 

for (10 s ) 


1.86 


0.620 

foL(10 s ) 


301 


316 

hi (A) 


1.10 


1.00 

hi (A) 


0.U4 


0.198 



0.6145 


0.6254 


manufacturer for each detector is a dose-equivalent-per- 
count, integrated over the energy spectrum of the calibration 
source. The calibration was transferred from a NIST M2 Cf 
standard using an Eberline “rem-ball” detector in conjunction 
with a Pu-Be source. This constant was used with 
measurements of the Wendi-II response function <%E) made 
by Olsher, et al ,, 42 to determine the response of the Swendi- 
II’s at the D-D neutron energy, 2.45 MeV. The uncertainty in 
the Swendi-Il calibrations is estimated, from data furnished 
by Olsher, et al ., and Eberline, to be 7.9%. For reference, S is 
one half of the D-D fusion reaction rate. 
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The Swendi-II detector is housed in a 23 cm 
diameter cylindrical moderator, 22 cm high, and both 
detectors were used tangent to the spherical chamber, in the 
equatorial plane (see Fig. 1). The observation angles 6 are 
measured from the transmitter viewport axis in Fig. 1, and 0, 
was fixed at 35° while 0 2 = 215° and 150°, respectively, for 
the 15 and 23 cm cathodes. Each detector is attached to a 
standard meter with a digital display, and the neutron counts 
S were recorded manually from these displays and logged 
into the data acquisition system where they are time-stamped. 
In some cases S was recorded directly, but for most of the 
data (*) the counts S were recorded and converted, post-test, 
to rates. At the measurement stations, the Swendi-II 
responses resulted approximately in one count per 1.4 x 10 4 
neutrons. Assuming a Jn stochastic signal-to-noise ratio, 
this gives a statistical error of 8%, equal to the calibration 
error, at 2.2 x 10 6 neutrons, which typically corresponds to an 
integration time of a few seconds. The thickness of the water 
jacket surrounding the vacuum chamber is 0.97 cm, and a 
simple calculation of the neutron emission signal attenuation 
due to this jacket gives a value on the order of a few percent. 
However, no correction for this effect is included in the data 
presented here. 

Neutron emission data from D-D fusion are given in 
Fig. 12. The data in Fig. 12(a) are compiled from several 
different tests overlaid on the time axis, and emission rates 
are plotted as functions of cathode voltage in Fig 12(b) for 
both cathodes. The cathode voltage was held approximately 
constant for the duration of each test in this series and the 
neutron rates are thus also approximately constant. Emission 
rates of a few MHz in the range 40 to 60 kV are typical of 
these devices. 2, 4> 25-27 The current was about 40 raA for the 1 5 
cm cathode at all voltages, but was reduced somewhat for the 
23 cm cathode, because of decreased stability of the plasma at 
high voltages for this cathode. However, normalization of the 
neutron rates to the current does not alter the qualitative 
features of Fig. 12(b). The dependence of S on cathode 
radius for these data is in contrast to that observed earlier at 
lower voltages in an “emitter-mode” 20 device, which showed 
an increase (about a factor of 2) in S between a = 4 and a - 
2. 4 The pressures in Fig. 12 are from = 0.52 to 0.4 Pa for both 
cathodes, compared to a fixed 0.25 Pa pressure in the emitter- 
mode device. 

As with the microwave data, comparison with the 
model is also of interest for the fusion data, although in this 
case the parameter space is less restricted, due in part to the 
absence of spatial resolution. Because of the high voltages, 
mono-energetic solutions with sufficient current are now 
possible, and for the dispersed distributions, only kinetic 
solutions are considered. Also, it should be noted that the 
visual appearance of the plasma at these voltages differs from 
that at lower voltages in that the spherical “core” is much 
smaller and in that distinct radial channels become 
discemable at grid openings. 22 

To calculate neutron rates from the model, values 
must be assigned to £, 77, T 2y and the atomic fractions Zo and 
Zi, defined as 



-F(kV) 

FIG 12. (a) Cumulative neutron counts for two directions for the 15 cm 
cathode Open symbols are #, solid symbols are (h (b) average emission 
rates as functions of cathode voltage, again for two directions, for the 15 cm 
cathode (upper curve) and the 23 cm cathode (lower curve). 


Xo 



Xi 


n D 


no # 

+"/>* * 


(29) 


where n p + , n p + , n D and are the number densities of, 

respectively, free deuterons, molecular deuterium ions, 
neutral deuterium atoms and neutral deuterium molecules. 
The secondary electron coefficient S depends on the ion 
energy, and thus on V g and fo, and from typical behavior of 
this coefficient 38 and earlier work, 39 a value of 8 = 1 is used 
here. For the currents 7 0 the efficiency is taken to be 17 = 1 . 
Both limiting values (0 and 1) of Xo are examined, and for 
convenience we set T 2 ~ 10 3 K and xi = 9. As indicated 
above T 2 is not known, and by the same token Xi may be 
substantial (certainly it is not categorically negligible, as 
evidenced by the Balmer-series optical emissions from the 
plasma), but the target deuteron number density, and thereby 
the fusion rate, is taken to be proportional to (2 - x 2 )T { 1 , and 
the results are thus readily generalizable with respect to T 2 
and X 2 * Beam-beam reactions are neglected. 4 Given these 
parameters, kinetic solutions for the experimental conditions 
of Fig. 12(b) were used to calculate neutron emission rates, 42 
and these rates are compared with the empirical rates from 
the figure, averaged over the two directions, in Table VI. The 

increase in S with increasing Xo is due to the increase in the 
fusion cross section <t f {E) with increasing incident deuteron 


energy. For the thermal solutions, e 0mim = -7, e l ^ mx = 8 was 


used for comparison; with U 0 ^ \q e V\ the Xo = 0 rates are 
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similar for the thermal and parabolic cases (*). These 
solutions are discussed further in the following section. 

Neutron emission was also observed, at reduced 
rates, in protium experiments conducted subsequent to 
operation with deuterium. This emission is attributed to 
residual deuterium from the vacuum chamber wails and 
cathode grid (*); the fusion cross section for protium at these 
energies is negligible. Fig. 13 shows emission rates for the 15 
cm cathode for four sets of measurements in protium, 
sequential in time, made on the 2nd, 3rd, 12th and 15th days, 
respectively, following termination of deuterium flow 
through the chamber. No operations were conducted on the 
intervening days. For comparison, background rates, 
measured with the high voltage switched off, are typically of 
order 1(T* MHz. 

Significant anisotropy in neutron emission was 
observed in the residual experiments. The anisotropy A is 
quantified here as the fractional difference between the two 
directions, 

A s 2 [S(4 ) - S(0 2 )]/[S(tf, ) + S(0 2 )] , (30) 

and is plotted in Fig. 14 for the individual measurements from 
which Fig. 13 was compiled, as well as for an analogous set 
of measurements made with the 23 cm cathode. The mean 
anisotropies were 45% and 16% for the 15 and 23 cm 
cathodes, respectively. By comparison, the measurement 
error for this data AS, computed as the root-sum-square of the 
calibration error and the statistical error at each point, has 
mean values of 12.6% and 8.7%, respectively. A similar 
analysis was performed for a collection of deuterium data 
(including that in Fig. 12) and these results are also shown in 
the figure. For normal deuterium operations, the anisotropy in 
neutron emission is within the measurement error. Since the 
cathode grid itself is somewhat anisotropic, both by design 
and construction, these results suggest that the contribution to 
S of beam-target reactions in the cathode is more significant 
under the residual conditions. The cathode contribution in 
normal operations, not addressed in Table VI, is discussed in 
Ref. 4. 

While the IEC device is generally characterized by 
stable operation, anomalies were observed on occasion. In 
one of the experiments, emission rates as much as a factor of 
2 greater than those of Fig. 12 were measured under nominal 
conditions approximately identical to those of the figure. 



during protium operation with the IS cm cathode. 



o.o 


0.1 0.2 
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FIG 14. Neutron emission anisotropy for the two cathodes, (a) Residual 
emission from protium operation. The mean value of |Aj /SS is 4.39 fix the 
15 cm cathode and 1.94 for die 23 cm cathode ( |A|/£S = 3.59 and 1.85). 
(b) Normal deuterium emissions. The mean value of |Aj /<ff is 0.562 for the 
15 cm cathode and 0.891 for the 23 cm cathode ( | A |/<W = 0.562 and 
0.896); i.e., the anisotropies are within the measurement error. 


Table VI. Measured and calculated D-D fusion neutron emission rates for 5= 77 = 1 The measured rates are denoted S , and the theoretical rates are denoted 
S E , S L and S T , respectively, for the mono-energetic, parabolic, and thermal distributions. The experimental conditions, and the measured rates, are those of 
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K^kV) 

P( Pa) 

A o(A) 

S 

4 

4 

S r 

4 

4 

4 


-38.3 

0.523 

0.573 

034 

0372 

0.109 

0.263 

5.24 

133 

233 

15 cm 

— 43. 1 

0.493 

0.545 

130 

0.585 

0.170 

0360 

630 

2.43 

2.74 

grid 

-50.5 

0.455 

0.572 

2.72 

1.09 

0333 

0.588 

10.2 

3.79 

3.90 


-61.4 

0.413 

0.573 

532 

2.12 

0.672 

0.981 

14.9 

5.91 

5.52 


—42.0 

0.517 

0.652 

039 

0.909 

0.272 

0.594 

11.1 

3.92 

433 

23 cm 

-49.4 

0.480 

0.616 

1.12 

1.66 

0.495 

0.894 

153 

5.66 

5.91 

grid 

-55.8 

0.400 

0.529 

1.17 

1.96 

0.588 

0.948 

15.1 

5.68 

539 


-62.9 

0.403 

0.452 

1.49 

3.18 

0307 

1.16 

173 

6.68 

6.23 
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Similarly, mean emission anisotropies of approximately 24% 
were noted in deuterium plasmas during one test sequence. 

VI. DISCUSSION AND CONCLUSIONS 

The electron density profiles measured in the IEC 
plasmas are generally consistent with dispersed total energy 
distributions which have thermalized angular momentum, i.e., 
distribution functions which vanish as the transverse kinetic 
energy goes to zero. The particular distributions used here are 
(1) that of thermal equilibrium, optionally truncated by step 
functions, and (2) a parabolic distribution function, 
computationally simpler than the thermal case and equivalent 
to high-temperature equilibrium truncated by the local 
potential. The magnitude of the mean densities for these 
distributions is limited by net charge near the cathode to 
values which are lower that those observed in most cases. For 
Ar at 1 kV, for example, the disparity is more than an order of 
magnitude. Solutions of sufficient density, however, are 
possible using the plasma approximation, and these solutions 
have transverse profiles which, like die Ar data but unlike the 
kinetic solutions, tend to increase near the cathode. 

Symmetric minima in the measured transverse 
profiles are observed at about r g 12 for some of the data. This 
feature is particulaiy clear for the standard 15 cm grid, a 
result which may, however, be due more to the transverse 
scan limits and the measurement uncertainties than to grid 
size or design. Similar but distinctly broader and shallower 
minima are predicted by 1-D solutions which superpose low- 
density quasi-mono-energetic kinetic distributions and high- 
density neutral thermal distributions. These hybrid solutions 
exhibit a double potential well in which the inner well has a 
depth of about 1/3 of the cathode potential. In all the 
experimental cases, though, it is worth emphasizing that the 
solutions are not generally determined uniquely by the data, 
and that where the agreement is limited, as in the hybrid 
case, a candidate cause is the limitations of the model; the f 
here are significantly restricted in form and are not explicitly 
coupled to coilisional processes in the plasma. 

In Figs. 10 and 11, theoretical transverse profiles 
from the 1-D model are matched to the empirical data. Based 
on the mean density of chained particles from these solutions 
and, again, with T 2 - 10 3 K and Xi = 0, the mean ionization 
fraction for the Ar plasmas is about 2 - 2.3 x lO^ 4 for either 
cathode and lower, from 0.3 - 1 x 10^ for the D 2 . These 
values, which include the hybrid case with its substantial 
radial structure, are veiy similar to those obtained from the 
mean measured areal densities assuming a uniform radial 
distribution. 

The effective grid transparency predicted by the 
high-density (plasma approximation) thermal solutions 
depends on the lower limit imposed on the total ion energy 
£o,mm (upper bound on £o). By reducing the slow, short- 
trajectory, ion populations near the grid, the current is 
reduced disproportionately to the total ion population, and 
absorption efficiencies near unity are possible. For 
unrestricted energies, both the thermal and parabolic 
distribution functions result in typical efficiencies of 0.1 to 


*.„< £ ,- |kev > 

0.25 0.5 0.75 1 



Fig. 15 Electron number density as a function of cathode voltage for the Ar 
pressure scans. The data are from Fig 5 and the number densities are 
calculated for uniform spatial distribution of the electrons. Is = 45 mA for 
all data. The thick curve shows the perveance scaling of density and the 
broken line, referenced to the top axis, illustrates the dependence of density 
on the lower bound of ion energy for fixed current in the thermal 
distributions. 


0.015. Exterior to the cathode, it is necessary to impose an 
upper bound < E g on the ion energy in order to obtain 
consistent solutions for the experimental conditions. 

In the Ar pressure data of Fig. 5(b), Ni(x=Q) 
increases with decreasing E g until about 0.75 keV, where 
jVi(£g) becomes approximately constant. Fig. 15 shows this 
explicitly, giving n x ( E g ) for the Ar data from both grids in 
Fig. 5 using a uniform radial profile; the current I s = 45 mA. 
Dispersion of the ion energy due to ionization distribution 20 
and charge exchange 44 has been used to explain a dependence 
of fusion rate density on neutral pressure resembling the n { 
curves of Fig 15 (cf. Fig. 4 of Ref. 20), and similar 
considerations apply here. From Eqs. (9) — ( 11), for a given 
solution </{p) and fixed current /o, the ion density scales as 
E g m and this is illustrated in the figure using an arbitrary 
= 4xl0 10 cm" 3 keV 3/2 . Assuming that n t = n Q9 which is 
generally the case for the kinetic solutions considered and 
necessarily the case for the high-density solutions, n, also 
will tend to increase with pressure like E g m . However, for 
the voltages of this data the charge exchange cross section a q 
for Ar + ions in Ar increases with decreasing energy, 
increasing by about a factor of 2 (to = 3 x HT 15 cm 2 ) between 
2.5 and 0.1 keV. 38 From 2.5 kV down to 0.75 kV the mean 
ion velocity U 0 decreases by about a factor of 2 while the 
neutral number density increases by about a factor of 10 
(see F ig. 5). Thus the charge exchange rate per ion n ^ 
o q u 0 will increase with decreasing voltage and Jo will tend to 
increase at the low energy end. Fig. 15 also gives a plot of 
n Q (£o^nax) for constant I 0 at E g = 1 keV and it is seen that 
n 0 decreases sharply as the low energy region is filled. The 
solutions used are those of the hybrid case from Fig. 1 1 with 
77 = 1 . Thus for distributions of this kind there are opposing 
influences on n x with increasing pressure and it is reasonable 
to expect that a local maximum like that in Fig. 15 would be 
observed. 
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For the parabolic distributions, on the other hand, 
variation of s> flTiax (and is largely ineffective in 

regulating I 0 independently of ^ . In this case, as in die 
thermal case with = U solutions consistent with the 
observations require variation of 7 with cathode voltage 
and, as indicated, imply small values for this parameter. 
These two alternatives, high-current/low-energy and low- 
current/high-energy, may, for a given mean density, be 
distinguished observationally on the basis of the interior 
density profile near the cathode. Both n and N in this region 
are larger for the high-current cases, as illustrated in Figs. 9 
and 10. 

For the deuterium fusion data, the measured quantity 
is the total neutron emission rate S . The measurement of 
phase shifts below 1 Pa in D 2 was not practical with the 
current system [e.g., see Fig. 5(c)], and consequently the 
densities are not directly known. In general though, the 
kinetic solutions for the small perveances of the 

combinations of 7 0 and V g for this data form shallow wells 
with V ~ V g , and 7? 0 in these cases is fairly insensitive to n x . 
The k x for the solutions of Table III were chosen to give /?, of 
the same order as n 0 , but the resulting values of S are 
similar in the absence of electrons. The presence of electrons 
in the solutions does, however, drive the central potential 
further towards the cathode potential. For the data in the 
table, 1 - (p = 0.001 to 0.005 is typical. These shallow wells 
are consistent with the argument given by Meeker, et al. y for 
low-energy electron populations trapped inside the cathode 
for ion-injected IEC plasmas (cf. Fig. 14 of Ref. 28). 

The low-energy electrons are a source of both 
atomic and molecular ions, in the range of from s 50 to 100 
eV, from H 2 data 45 and from vibrational state tables for D* 46 
the cross section for electron impact ionization of D2 is = 
10 -16 cm 2 . Also from the tables, the cross section for 
dissociative ionization of D2 (D2 + e -> D -1- D + + 2e) may be 
estimated as also ^ 10“ 16 cm 2 . 46 In these conditions, then, it 
is expected that both n /r and n D + may be appreciable, and 
from Table VI % 0 >0 is consistent with the S measurements 
at higher voltages for the 15 cm grid. For these data, though, 
the calculated values of S increase with the grid radius, as 
shown in the table, and this was not observed. It may also be 
noted that in the present case ionization by ions is more 
significant than in Meeker's case, due to the larger cross 
section. In the range 30 to 60 kV, the cross section for proton 
impact ionization of H 2 is =2 xl0~ 16 cm 2 , similar to the 
charge exchange cross section. 47 

The spatial distribution of the fusion rate density 
F(p) in the dispersed solutions is roughly constant inside the 
cathode and decreases quickly outside the cathode 
(independently of XoY The mono-energetic solutions in 
contrast have pronounced central peaks. However, all three 
forms in Table VI have similar cumulative profiles 
C(p ) « F(p' } )p' 2 dp' and inside the cathode these 

profiles resemble that given in Fig. 4 of Ref. 4 from proton 
collimation measurements. For radii well outside the cathode, 
on the other hand, the empirical profile increases 
monotonically while the kinetic solutions give flat due 
to w o ,^-»0 . Proton collimation results reported in Ref. 22 
show a pronounced central peak in F(p ) . Also, it may be 


remarked that, for fixed J 0 in the thermal distributions, the 
behavior of n Q (£ 0 ^) suggests that any decrease in with 
increasing E g which is less marked than the scaling may 
be of fusion interest, i.e., due to increase in U 0 /Eg. 
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